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ABSTRACT 

It is shown that detecting or setting an upper hniit on the scalar gravitational 
radiation is a good experimental test of relativistic gravity theories. The relativistic 
tensor-field theory of gravitation is revised and it is demonstrated that the scalar 
monopole gravitational radiation must be added to the usual quadrupole radiation. 
In the case of the binary pulsar PSR 1913 + 16 it is predicted the existence 0.735% 
excess of the gravitational radiation due to the scalar gravitational waves. 

1. Scalar Gravitational Radiation as a Test of Gravity Theories 
1.1. Geometry of Spacetime or Quantum Field in Flat Spacetime 

Gravitational wave experiments give us a new possibility in choosing among 
alternative gravity theories. Here we revise nonmetric tensor-field theory (TFT) 
of gravitation, i.e. relativistic theory of symmetrical second rank field ip^'^ in fiat 
Minkowski spacetime 77*^ and compare its gravitational wave predictions with general 
relativity (GR) ones. 

It is easy to show that widespread opinion on full coincidence of TFT and 
GR (see for example^) based on the assumption of the uniqueness of the energy- 
momentum tensor (EMT) of the gravitational field, which was used for the iteration 
procedure. However, in the framework of Lagrangian formalism of relativistic field 
theory, EMT of any field is not defined uniquely: T*^ =^ T*'' + ^^'"'^ for = -<|)*'^ 
(see for example^' ^), and one needs additional physical restrictions to choose final 
form of the EMT (for example such conditions as energy positiveness, tracelessness, 
symmetry). It is apparently that different EMTs will lead to different nonlinear 
theories of gravitation. 

To illustrate the above-said discussion let us consider nonlinear generaliza- 
tion of Poisson's equation for the case of distributed source describing negative and 
positive energy density of the gravitational field 



and 



The solution of Eq.la is 



= -(V(^)Vc^ (la) 
Acp = +(V(^)Vc'. (16) 



~ V ' " r2(l-GM/c2r)' ^ ^ 

whereas the solution of Eq.lb is 

^Prom Eq.2 and Eq.3 we see two possible ways for the construction of the nonlin- 
ear TFT. The first one, based on the negative energy density of the gravitational 
field, leads to the infinite gravity force at the finite distance Rg = GM/c^. The 
second way, based on positive one, gives TFT without singularity. 

As for GR, the energy density of the gravitational field is a poorly defined 
concept as a consequence of geometrical interpretation of gravity (there is no tensor 
characteristic of the gravitational field EMT). For example, Landau-Lifshiz pseu- 
dotensor gives = — 7(V(/9^)^/87rG' but Grischuk-Pctrov-Popova one gives t^^ = 
— llCVipj^y/STrG for the spherically symmetric static (SSS) weak field in harmonic 
coordinates. 

Within the scope of TFT there is a real tensor quantity for the energy density 
of the gravitational field. For example canonical EMT gives T^^ = +(Vv9^)^/87rG 
for SSS weak field and corresponds quantum description of the gravitational field as 
aggregate of gravitons in fiat spacetime. Gravitons are massless particles, i.e. some 
kind of matter in spacetime, which carry-over positive energy and momentum in the 
spacetime. Besides the sum of two tensors ■0*'^ + V^'^ — 9^^ is not a metric tensor 
because the covariant components -01^ -|- r]ik — Qik of this tensor provide the mixed 
components ip], + r]l = gl ^ Si and the trace g^^Qik = 4 + 2-?/' + 0{il)ff}) ^ 4. It means 
that TFT is a scalar-tensor theory, but not pure tensor one, and hence includes spin 2 
and spin gravivtons. 



1.2. Post-Newtonian Tensor-Field Theory of Gravitation 

Let us consider the relativistic symmetric tensor field ■0**^ in Minkowski's space- 
time rf^. As so as really observed gravitational fields are the weak ones (|(/?| <^ c^), it 
is naturally to begin the construction of TFT for the weak field case. In this case we 
have a very closed analogy with the electromagnetic field and can use the standard 
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Lagrangian formalism of the relativistic field theory (bellow we utilize notations of 
the text-book^). 

We begin with the action integral in the form 

^ ^ ^{9) + ^(int) + S(p) = - J{A(^g) + A(^int) + A(p))(iQ (4) 

where (g), {int) and (p) indicate gravitational field, interaction and particles parts of 
actions and Lagrangians. The Lagrangians are given by the following expressions 

A 1 



Mint) = -\i^lmT[^), (6) 



= -V^kTl^y (7) 

It has been shown by Kalman^ and Thirring^, that the total EMT of the 
system contains three parts, which correspond to ones of the action integral (4), 

^ (E) - ^ (p) + ^ {int) + ^ {a) \P) 

where the canonical EMT of the gravitational field for Lagrangian (5) has the form 



the interaction EMT is 



- l^J^'" - ^jS)^lmU'u-, (10) 



the point particles EMT is 



Tip) =T.mac'S{r - r„){l - -|}^/X«a- (H) 



Therefore the nonlinear TFT must include the interaction Lagrangian in the form 

1 

c 



hint) = -^^ImTl^)- (12) 



But the weak field condition allows us to use linear approximation as the first step 
and then to make nonlinear corrections (Post-Newtonian TFT). 

The variation of the gravitational potentials in the action integral (4), where 
for fixed sources in the PN approximation we can use the interaction Lagrangian (12), 
yield the PN field equations in the form 

^ ^u,k ^ ^ki, _ ^,ik _ rj^k^im^^^ + ^ik^,i ^ ^ Tik^ (13) 
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Eq.l3 automatically requires conservation of the total EMT and leads to the 
particles motion equations in the form T^^^ fe = 0- 

The field equations are invariant (for fixed sources) under the gauge trans- 
formation ijj^'^ =^ ijj^^ + + 6^'^ and one can achieve Hilbert gauge in the form: 
'^^'^ k = iV''*- ^^i^ fi^^'i Eq.l3 become 



nijj 



ik 



SttG 



rpik 



(14) 



The very important feature of Eq.l3 is the multi component structure of the 
gravitational field description. Initial symmetric tensor ip^'^ simultaneously describes 
a mixture of four particles 



{t/ji'^} = {2} e {1} e {0} e {o'} 



(15) 



with spins equal 2, 1, and . It corresponds to 10 independent components of the 
tensor ■0*'^. After four gauge conditions, Eq.l4 describes a mixture of two particles 
which have two sources 



{^■'^■} = {2} © {0} ^ {r^} = {2} © {0} 
So that eq.(14) can be written in the form 



SttG 



ik 
{2} 



or ad,''' = — — 



ik 



rp 

4^ Hp) 



(16) 



(17a) 



and 



SttG 



ik 

2 - {0} 

l^ik„ 



or 



SttG 



ik 



(176) 

besides 



4 c2 "''A 

where ip^'' = ip)^2} +'^{0} ~ ^^'^ + ^if^'ip (and the same representation for Ti^y,))^ 
we take into account that gravitons of both kinds are massless particles, i.e. have 
traceless EMTs (T(s) = Tp). 

The variation of action integral (4) with respect to particle coordinates gives 
the equation of motion in the gravitational field^' ^ 



where 



Al = 



ds 



-Biiu^u^ 



2 2 



Bil = -^iPli- -^iPkt - -^lPkl,nU''u\ 



(18) 

(19) 
(20) 
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1.3. Static Weak Field 



In the case of spherically symmetric static (SSS) weak field the first approxi- 
mation EMT has a very simple form T^*|-) = diag(poc^, 0, 0, 0) and solution of Eq.l4 
is Birkhoff 's potential 

V'^'= = (^,diag(l, 1, 1, 1), (21) 

where ^p-^ is Newtonian potential (</7j^ = —GM/r outside the gravitating body). 

The gravitational field (21) can be expressed as the sum of the spin 2 and 
spin parts 

V'^'= = V'g}+^fo} = (/.,diag(3, 1, 1, l) + (-(^Jdiag(l, -1, -1, -1). (22) 

Using the SSS solution (21) and corresponding EMT expressions (9), (10), (11) 
we find the total energy densities of the system from Eq.8 

= T^) + + = (Poc^ + e) + po^. + (V^J^ (23) 

where (poc^ + e) is the rest mass and the kinetic energy densities, po<^N is the inter- 
action energy density, V^^^/SttG is the energy density of the gravitational field. The 
total energy of the system will be 

^(s) = / T^i)dV = Eo + Ek + Ep (24) 

where Eq — J{poC^) dV is the rest-mass energy, E^ — /(e) dV is the kinetic energy, 
and Ep is the classical potential energy that equals the sum of the interaction and 
the gravitational field energy: 

Ep = E^,^t) + E^,) = J {poif^ + -^(Vif^y) dV = ^l poif, dV (25) 

It is important to note that the canonical EMT (9) include both EMTs for 
spin 2 and spin gravitons. Hence the T^g^ given by (9) is the sum of energy densities 
of both components (see for details^). 

For external field of SSS body in PN approximation the only needed correction 
is to take into account the energy density of the gravitational field. Hence we have 
exactly Eq.lb and nonlinear addition for Birkhoff 's component is 

r = + (26) 

Substituting Birkhoff 's potential (21) into (18) and taking into account non- 
linear PN correction (26) one gets the three dimensional equation of motion for test 
particle 

-._(,4,4^)V,„.4I(I.V.„). (27) 
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For Eq.27 the pericenter shift of test particle is 



in which the nonhnear contribution (26) provides 16.7% of the value (28). Therefore 
in TFT there is a direct testing of the positive energy density of the gravitational 
field. 

It is easy to show^' ^' ^ that all PN classical relativistic gravitational effects have 
the same values as in GR but the other interpretations. For example substituting 
Eq.22 in Eq.l8 one finds that usual Newtonian force is the sum of the attracting force 
(spin 2) and repulsing force (spin 0) 

F = F{2} + F{o} = ~moV(f^ + ^moVifi^ = -moVifi^ = (29) 



1.4- Free Field 

Within the scope of TFT not only static, but also variable gravitational field 
contains the sum of spin 2 and spin gravitons. The field equations in the form (17), 
can be written for the free field as 

□(/.^'^-O, 0% = O (29a) 

and 

^n^jjT)''' = 0. (296) 
Eqs.29 can be derived from the free field Lagrangians 

A{2} = Y^0im,n0'"^'", (30a) 

and 

Eqs.30 is a consequence of the gravitational field Lagrangian (5) for the field connected 
with sources (full field equations (13)). Note that sign of (30b) is positive due to 
positive energy density condition for spin free particles. Corresponding EMTs for 
tensor (spin 2) and scalar (spin 0) gravitational waves are 

rg} = 

and 
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2. Emission of Tensor and Scalar Gravitational Waves 



Let us consider the retarded potentials solution of Eqs.l7 in the wave zone for 
slow motions in the source. For Eq.l7a we have usual quadrupole radiation (spin 2 
gravitons) with luminosity 

L{2} = (32) 

where Dap is the reduced quadrupole moment of the system of particles. 
In the same conditions the solution of the Eq.l7b is 

, 2GMo 2GEk 2GMo ^ G r 

V^(r, t) ^ \ In • R) + -^HaUpIap, (33) 

where Mq = ^vHa is the rest mass, E^^ = l/2J2fnaVa is the kinetic energy, R = 
J2 ^a^a/ S i^a IS the Center of mass, I^p = f^a^'^x^ is the moment of inertia of the 
system. Differentiation with respect to the time of Eq.33 leaves the only monopole 
term 

^(r, t) ^ (34) 

because the first term equals to zero in consequence of rest mass conservation, the 
third term equals to zero in view of inertial motion of the center of mass and the fourth 
term vanishes so that it has been taken into account within quadrupole radiation (32). 
Substituting Eq.34 into scalar wave EMT (31b) and averaging it on the sphere one 
gets the scalar monopole luminosity (spin gravitons) 

^{0} = ^,El (35) 



3. Detection of Tensor and Scalar Gravitational Waves 

3.1. Test Particles in Tensor Wave 

Let us consider inertial frame of reference where we analyzed action integral 
(4) and derived equations of motion for test particles (18). We can fix Cartesian 
coordinates with being the gravitational radiation source in the origin and the x-axis 
along propagation of gravitational wave. For plane tensor wave (jf^ satisfying Eq.29a 
the only nonzero components are = —0^^ = and (jP'^ = = Ax, where for 
monochromatic wave A{t, x) = cos{ujt — kx). Substituting above expressions into 
Eq. 18 and leaving the main terms one gets equation of motion of test particle in the 
plane monochromatic tensor wave 



It follows from Eq.36 that the tensor wave is transversal, i.e. there is no acceleration 
along propagation direction. Besides for the rest test particle {v^ — Vy = = 0) 
there is no interaction between particle and tensor wave (in a manner like the rest 
charged particle has no interaction with magnetic field in electrodynamics). Therefore 
in the case of tensor gravitational wave the velocity of test particle plays the role 
of "gravitational charge" and we need such detectors in which there is system of 
particles with nonzero relative velocities. For example it may be rotating bodies and 
even metallic bars where the degenerate electrons have nonzero velocities relatively 
crystallic proton lattice. The amplitude of the velocity deviation can be estimated 
as A^; = 2vA/c^. 



3.2. Test Particles in Scalar Wave 

For the scalar plane monochromatic gravitational wave satisfying Eq.29b we 

have t/'iqi = yl(t, x)ini'^^ = cos{ujt — kx)r]^'^. If we substitute this expression into 
Eq.l8 and leave the main terms we get the following equation of motion of test 
particle in scalar wave 



dvx A 
11 = '-- 



^ = 
dt 

(37) 

dt 

According to Eq.37 the scalar wave is longitudinal and the rest test particle 
oscillates around initial position with the velocity amplitude Ai> = cAjc? and the 
distance amphtude Ax = A/k(?. For two particles at a distance Iq <^ \ along x-axis 
we get the dimensionless amplitude of oscillation in the form AIq/Iq = A/c^. Hence 
the detection of the scalar wave can be achieved by means of bar detectors and 
laser interferometric detectors. It is important to note that scalar gravitational wave 
does not interact with the electromagnetic field because the interaction Lagrangian 
is ^g)r,, = \^l^VikTl^^^ = 0. 



4. Binary Pulsar PSR 1913+16 

4.1. Gravitational Radiation from Binary System 

For a binary system the loss of energy due to tensor gravitational radiation is 
given by Eq.32. Direct calculations lead to the well-known formula for quadrupole 
luminosity of the binary system 

< E >{2}= , ^ (38) 
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where mi, 1712 are masses of the two bodies, a is the semimajor axis and e is eccen- 
tricity of the relative orbit. The similar calculations based on Eq.35 give the loss of 
energy due to scalar monopole radiation^ 

G^m?mn (mi + m.2) ( + 76"^ ) 
< E >{o}= \ (39) 

The ratio of the scalar (39) to the tensor (38) luminosity is 

< E >{o} _ 5 (e' + i^^) 



<E>{2} 128 (l + Z|e2 + |e^ 



(40) 



Note that the value of the ratio (40) lies in the interval [0, 1.1%] and for the circular 
orbit equals to zero. However for a pulsating spherically symmetric body there is no 
quadrupole radiation and scalar radiation becomes decisive. In particular it follows 
from this that in the TFT it is impossible to have a "quiet" relativistic collapse of a 
spherical body. 



4.2. PSR 1913 + 16 

According to^'^° for the binary pulsar PSR 1913 + 16 the eccentricity e = 
0.6171309(6), hence the expected scalar radiation part (40) is 0.735 %. As far as the 
rate of orbital period change P is proportional to the total energy loss (< P >= 
— 1^ < E >), one expectes the corresponding excess of orbital period decrease due 
to scalar gravitational radiation. The observed value of the orbital period change is^O 

p(obs) _ _2.425(10) • 10-^^ 
while the theoretical prediction for pure tensor quadrupole radiation is^ 

p|thcor) ^ _2.402576(69) ■ lO'^^. 
Therefore the observed excess of the orbital period decrease is 

((Obs) - (Theor)) = +0.96 % ± 0.4 % 

that is very close to the expected value 0.735 % for the scalar gravitational radiation. 
Taking into account the spin part (39) of the gravitational radiation we get the 
theoretical prediction 

p(theor) _ p(theor) ^ p^theor) _ _2.420254(69) • 10"^^ 

in a good agreement with the observed value. 
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4-3. Problem of Galactic Acceleration 

It has been shown in^ that one must take into account the effect of the galactic 
acceleration of the pulsar and the Sun, and that of the proper motion of the pulsar. 
Note that the distance d to PSR 1913 + 16 is the very sensitive parameter for the 
calculation of the galactic effect. Unfortunately the line of sight to the pulsar pass 
through very complex region of our Galaxy and one must be very careful to use 
known distances to other pulsars for distance estimation to PSR 1913 + 16. Indeed 
in^ it was used indirect arguments to re-estimate the standard dispersion-measure 
distance of 5.2 kpc and to get the new distance of 8.3 kpc. For the d — 8.3 kpc the 
galactic effect is -1-0.69 %. However another arguments^^ based on an analysis of the 
pulse structure of PSR 1913-1-16 itself lead to an estimation of the distance about 
3 kpc. For the d = 3 kpc the galactic effect is +0.11 %. In any case it is clear that 
we need further investigation of the pulsar and the direct estimation of the distance 
to PSR 1913 + 16. 
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